Solidity of viscous liquids II: Anisotropic flow events 
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Recent findings on the displacements in the surroundings 
of isotropic flow events in viscous liquids [Phys. Rev. E, in 
press] are generalized to the anisotropic case. Also, it is shown 
that a flow event is characterized by a dimensionless number 
reflecting the degree of anisotropy. 
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In a previous paper henceforth referred to as (I) jjj , it 
was argued that viscous liquids close to the glass transi- 
tion P - where viscosity is roughly 10 15 times larger 
than that of, e.g., room temperature water - are more like 
solids than like the less viscous liquids studied in stan- 
dard liquid theory Jl2|,^3|. The idea || that viscous liq- 
uids are qualitatively different from less-viscous liquids is, 
of course, not new. It is a rather obvious idea, given the 
following fact. While "ordinary" less- viscous liquids have 
relaxation times in the picosecond range, i.e., comparable 
to typical phonon times, viscous liquids have much longer 
average relaxation times (roughly given by Maxwell's ex- 
pression r = rj/Goo, where r\ is the viscosity and G M the 
instantaneous shear modulus). This decoupling of relax- 
ation times from phonon times is also reflected in a decou- 
pling of diffusion constants |Q: For less- viscous liquids 
the molecular diffusion constant D is of the same order 
of magnitude as the transverse momentum diffusion con- 
stant, the dynamic viscosity of Navier-Stokes equation 
v oc rj. However, with increasing viscosity D decreases 
(roughly as r/^ 1 from a simple Stokes-Einstein type ar- 
gument) while v increases. At the glass transition v is 
about I0 30 times larger than D. 

The average relaxation time increases dramatically 
upon cooling. Goldstein has argued jl5| that already 
when r becomes longer than about f nanosecond is there 
a gradual onset of typical viscous-liquid behavior. As 
noted first by Angell §, this is roughly at the tem- 
perature below which ideal mode-coupling theory fllf | 
breaks down. It is generally believed that in viscous 
liquids "real" molecular motion beyond pure vibration 
takes place on the time scale defined by r, although in- 
homogeneities are likely to give rise to faster relaxations 
in some parts of the liquid ff7|-|22|| . "Real" motion is rare 
because it involves overcoming energy barriers large com- 
pared to fc^T J|Jf^]- The transition itself is a jump be- 
tween two potential energy minima, a process that lasts 
just a few picoseconds. One thus arrives at the follow- 
ing picture: Most molecular motion in a viscous liquid 
is purely vibrational, real motion is rare and takes place 



via sudden molecular rearrangements. It is interesting 
to note that this old picture J2|,^3| has never really been 
challenged (while the nature of the energy barrier to be 
overcome in the transition is still being debated p4|). 
In fact, extensive computer simulations have now defini- 
tively confirmed the picture pq ] . 

The sudden molecular rearrangements [p|]26| |3l| ] are re- 
ferred to as "flow events" below. It is generally believed 
that flow events are localized in the sense that only a 
limited number of molecules experience large displace- 
ments, while all other molecules are only slightly dis- 
placed; the large-displacement molecules involved in a 
flow event thus define a "region" of the liquid. Because 
flow events are rare and molecules most of the time just 
vibrate, a viscous liquid looks much like a solid. In (I) 
the small displacements in the surroundings of a flow 
event were calculated from solid elasticity theory assum- 
ing spherical symmetry. It was shown that the displace- 
ment u in the surroundings of a region is given by (where 
r is the distances to the region) 

u oc — . (I) 

The displacement is purely radial. However, spherical 
symmetry of flow events is not realistic; when molecules 
move from one potential energy minimum to another 
there must be some violation of spherical symmetry, even 
if the molecules were spheres with only radially depen- 
dent interactions. One is thus lead to ask to whether Eq. 
(|l|) and its consequences remain valid in the anisotropic 
case. 

As in (I) the starting point is the solidity of viscous 
liquids as reflected in the slow "real" motion of the 
molecules. This fact implies that the average force on 
any molecule is extremely close to zero. In a continuum 
description, the average force per unit volume is the di- 
vergence of the stress tensor cry, where i,j — 1,2,3 are 
spatial indices. The condition of average zero force - elas- 
tic equilibrium - is (where di = d/dxi and one sums over 
repeated indices) 

diCTij = 0. (2) 

Linear elasticity theory |}2) may be applied to the re- 
gion surroundings, because the molecular displacements 
in these surroundings are small and because there is elas- 
tic equilibrium in the liquid before as well as after a flow 
event. Most likely, there are large "frozen-in" stresses in 
the liquid, but the changes in the stress tensor induced 
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by one flow event are small, except in the region itself. 
Now, define a sphere centered at the region, large enough 
that outside the sphere the flow event induced displace- 
ments and stress tensor changes are so small that lin- 
ear elasticity theory applies for the changes. Imagine all 
molecules within the sphere being removed and the forces 
from these molecules acting on the molecules outside the 
sphere replaced by external forces applied to the surface 
of the sphere. This is done before as well as after the 
flow event. The flow event induced displacements of the 
surroundings can then be calculated from the change of 
these external forces. To do this we first consider the 
distance dependence of displacements in an elastic solid 
when an external force is applied to just one point. There 
is then a continuous flow of momentum into the solid at 
that point. The stress tensor is the momentum current 
and the mechanical equilibrium condition Eq. (||) is the 
zero-divergence equation reflecting momentum conserva- 
tion. By considering Gauss surfaces at various distances 
from the point, one concludes from Eq. (||) that the stress 
tensor decays as r~ 2 , where r is the distance to the point. 
Since the stress tensor is formed from first order spa- 
tial derivatives of the displacement u, we conclude that 
u oc r -1 Q. This result applies also when several exter- 
nal forces are applied to the solid, as long as these forces 
do not sum to zero. In our case, the external forces re- 
placing the forces from the molecules within the sphere 
do sum to zero: The forces from the molecules outside 
the sphere on those inside must sum to zero - otherwise 
the latter molecules would start to move. By Newton's 
third law, the sum of the forces acting from the molecules 
inside the sphere on those outside - the forces that are 
replaced by external forces - must therefore also sum to 
zero. When the external forces sum to zero, the stress 
tensor does not decay as r~ 2 but as r~ 3 (the mathemat- 
ics behind this fact is the same as that implying that 
the electric field from a charge distribution with zero to- 
tal charge decays as r~ 3 and not as r~ 2 ). Consequently, 
since the stress tensor is given as first order derivatives 
of the displacement vector we arrive at Eq. ([j]), which 
is now to be understood as valid for each of the three 
components of the displacement vector. In particular, 
we note that the predictions of (I) for the displacement 
and rotation angle distributions in the surroundings of 
a flow event (P(u) oc uT^I 2 and P(4>) oc 4>~ 2 ) are valid 
also in the anisotropic case. The first prediction has re- 
cently been confirmed in computer simulations of a bi- 
nary Lennard- Jones mixture (34| , the second is consistent 
with the small rotation angle distribution tentatively in- 
ferred from NMR experiments by Bohmer and Hinze on 
glycerol, P(<j>) oc l/sin 2 (» Q. 

We now show that it is possible to characterize flow 
events according to their anisotropy. The elastic equilib- 
rium in the surroundings of a flow event region before as 
well as after the flow event implies that the stress tensor 
change, Aaij, has zero divergence (i.e., obeys Eq. (ph). 



Since u and Aaij are linked by linear elasticity theory 
one has |52| 

V 2 (V • u) = . (3) 

This equation can be solved asymptotically for r — > oo: 
Equation (Q) implies V • u oc r~ 3 . Any real solution to 
the Laplace equation decaying as r -3 can be written ]36| ] 
as a P2(0, 0)/r 3 , where a > is a constant and Pi is a 
normalized linear combination of second order spherical 
harmonics: P 2 = J2lnZ-2 c ™ Y i™' where c* n = c_ TO and 
J2m= 2 -2 l c «i| 2 = 1- This expression applies far away from 
the flow event: tq <C r, where r$ is the region size. On 
the other hand, the expression does not apply beyond 
the "solidity length" I discussed in (I), where essentially 
no flow event induced displacements are expected. Since 
V • u is dimensionless a has dimension (length) 3 . Writ mg 
a = a/po, where po is the average (number) density and 
a is dimensionless, we have 

V • u = a ( r „ « r « I) . ( 4 ) 

Po r J 

The parameter a is a measure of the flow event 
anisotropy, the case a = corresponding to isotropic flow 
events. 

In a homogeneous system described by linear elasticity 
theory the density change following an elastic displace- 
ment is equal to — poV • u |52|. Thus, if a viscous liquid 
were homogeneous the density change in the surround- 
ings of a flow event would be given by Eq. (|J) (looking 
like an electronic d-orbital). However, the density of a 
viscous liquid is not quite spatially constant. As is easy 
to show, the density change induced by a flow event has 
an extra term, — u ■ Vp, coming from the fact that the 
whole density profile is displaced. Far away from the 
flow event this extra density change term dominates over 
the (V • u)-term. 

The flow event induced changes given by Eqs. (|l|) 
and (Q) were calculated from the fact that there is a 
linear relation between displacement and stress tensor 
change. Therefore, these results are valid independent 
of the chemical nature of the liquid. One possible ob- 
jection to these results is that dynamic inhomogeneities 
most likely give rise to spatially varying elastic constants. 
However, being mainly interested in the high viscosity 
limit where the solidity length is large, these inhomo- 
geneities are not expected to have any significant effect 
on the average displacements in the surroundings of a 
flow event (the "long wave length" limit). Finally, we 
note that the sharp distinction between "real" motion 
and vibration is somewhat blurred by the fact that "real" 
motion takes place not only in the region itself in the 
form of large jumps but also in the surroundings in the 
form of small jumps. However, as may be shown from 
Eq. (pi), the dominant contribution to the mean-square 
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displacement of a molecule comes from the "real" motion 
of molecules inside regions. 

To summarize, arguing from the "solidity" of viscous 
liquids, the flow induced displacements in the surround- 
ings have been calculated for the general, anisotropic 
case. It has been shown that the r-dependence of these 
displacements is the same as that induced by isotropic 
flow events. A dimensionless number a has been intro- 
duced as a measure of the degree of anisotropy of a flow 
event. 
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